Reading in MappingMuseumData:


museum_df <- read_csv("MappingMuseumsData2021_09_30.csv")
Rows: 4191 Columns: 35── Column specification ────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (24): museum_id, Name_of_museum, Address_line_1, Address_line_2, Village,_Town_or_City, Postcode, Admin_area, Accreditation, Governa...
dbl (11): Latitude, Longitude, DOMUS_identifier, Area_Deprivation_index, Area_Deprivation_index_crime, Area_Deprivation_index_education,...
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
tidy_museum_df <- museum_df %>%
  pivot_longer(cols = c(starts_with("Area_Deprivation_Index")), names_to = "Area_Deprivation_Type", values_to = "Area_Deprivation_Score") %>%
  select(-ends_with("_provenance"),-starts_with("DOMUS"),-ends_with("_code")) %>%
  mutate(status = "open") %>%
  mutate(status = ifelse(Year_closed != "9999:9999","closed",status))
tidy_museum_sum <- tidy_museum_df %>%
  filter(Area_Deprivation_Type != "Area_Deprivation_index") %>%
  separate(col = "Area_Deprivation_Type", into = c(NA,NA,NA,"Area_Deprivation_Type"),sep = "_") %>%
  group_by(Area_Deprivation_Type,Area_Deprivation_Score,status) %>%
  summarize(n = n())
`summarise()` has grouped output by 'Area_Deprivation_Type', 'Area_Deprivation_Score'. You can override using the `.groups` argument.
  
ggp <- ggplot(tidy_museum_sum,aes(x = Area_Deprivation_Score,y = n, fill = status)) +
  geom_col(position = "dodge") +
  facet_wrap("Area_Deprivation_Type",) +
  scale_x_continuous(breaks = scales::pretty_breaks(n = 10)) +
  theme_classic()

ggplotly(ggp,tooltip = c("x","y","status"))
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0NCg0KbGlicmFyeSh0aWR5cikNCmxpYnJhcnkocmVhZHIpDQpsaWJyYXJ5KGdncGxvdDIpDQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeShzdHJpbmdyKQ0KbGlicmFyeShwbG90bHkpDQoNCmBgYA0KDQpSZWFkaW5nIGluIE1hcHBpbmdNdXNldW1EYXRhOg0KDQpgYGB7ciByZWFkLWluLWRhdGF9DQoNCm11c2V1bV9kZiA8LSByZWFkX2NzdigiTWFwcGluZ011c2V1bXNEYXRhMjAyMV8wOV8zMC5jc3YiKQ0KDQpgYGANCmBgYHtyIHRpZHktdXAtZGF0YX0NCnRpZHlfbXVzZXVtX2RmIDwtIG11c2V1bV9kZiAlPiUNCiAgcGl2b3RfbG9uZ2VyKGNvbHMgPSBjKHN0YXJ0c193aXRoKCJBcmVhX0RlcHJpdmF0aW9uX0luZGV4IikpLCBuYW1lc190byA9ICJBcmVhX0RlcHJpdmF0aW9uX1R5cGUiLCB2YWx1ZXNfdG8gPSAiQXJlYV9EZXByaXZhdGlvbl9TY29yZSIpICU+JQ0KICBzZWxlY3QoLWVuZHNfd2l0aCgiX3Byb3ZlbmFuY2UiKSwtc3RhcnRzX3dpdGgoIkRPTVVTIiksLWVuZHNfd2l0aCgiX2NvZGUiKSkgJT4lDQogIG11dGF0ZShzdGF0dXMgPSAib3BlbiIpICU+JQ0KICBtdXRhdGUoc3RhdHVzID0gaWZlbHNlKFllYXJfY2xvc2VkICE9ICI5OTk5Ojk5OTkiLCJjbG9zZWQiLHN0YXR1cykpDQpgYGANCg0KYGBge3Igdml6LXN1bS1vcGVuLXZzLWRlcHJpdmF0aW9ufQ0KdGlkeV9tdXNldW1fc3VtIDwtIHRpZHlfbXVzZXVtX2RmICU+JQ0KICBmaWx0ZXIoQXJlYV9EZXByaXZhdGlvbl9UeXBlICE9ICJBcmVhX0RlcHJpdmF0aW9uX2luZGV4IikgJT4lDQogIHNlcGFyYXRlKGNvbCA9ICJBcmVhX0RlcHJpdmF0aW9uX1R5cGUiLCBpbnRvID0gYyhOQSxOQSxOQSwiQXJlYV9EZXByaXZhdGlvbl9UeXBlIiksc2VwID0gIl8iKSAlPiUNCiAgZ3JvdXBfYnkoQXJlYV9EZXByaXZhdGlvbl9UeXBlLEFyZWFfRGVwcml2YXRpb25fU2NvcmUsc3RhdHVzKSAlPiUNCiAgc3VtbWFyaXplKG4gPSBuKCkpDQogIA0KZ2dwIDwtIGdncGxvdCh0aWR5X211c2V1bV9zdW0sYWVzKHggPSBBcmVhX0RlcHJpdmF0aW9uX1Njb3JlLHkgPSBuLCBmaWxsID0gc3RhdHVzKSkgKw0KICBnZW9tX2NvbChwb3NpdGlvbiA9ICJkb2RnZSIpICsNCiAgZmFjZXRfd3JhcCgiQXJlYV9EZXByaXZhdGlvbl9UeXBlIiwpICsNCiAgc2NhbGVfeF9jb250aW51b3VzKGJyZWFrcyA9IHNjYWxlczo6cHJldHR5X2JyZWFrcyhuID0gMTApKSArDQogIHRoZW1lX2NsYXNzaWMoKQ0KDQpnZ3Bsb3RseShnZ3AsdG9vbHRpcCA9IGMoIngiLCJ5Iiwic3RhdHVzIikpDQpgYGANCg==